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ABSTRACT 

We present radial velocity observations of four extremely low-mass (0.2 M & ) 
white dwarfs. All four stars show peak-to-peak radial velocity variations of 540 
- 710 km s _1 with 1.0 - 5.9 hr periods. The optical photometry rules out main- 
sequence companions. In addition, no milli-second pulsar companions are de- 
tected in radio observations. Thus the invisible companions are most likely white 
dwarfs. Two of the systems are the shortest period binary white dwarfs yet dis- 
covered. Due to the loss of angular momentum through gravitational radiation, 
three of the systems will merge within 500 Myr. The remaining system will merge 
within a Hubble time. The mass functions for three of the systems imply com- 
panions more massive than 0.46 M ; thus those are carbon/oxygen core white 
dwarfs. The unknown inclination angles prohibit a definitive conclusion about 
the future of these systems. However, the chance of a supernova la event is only 
1% to 5%. These systems are likely to form single R Coronae Borealis stars, 
providing evidence for a white dwarf + white dwarf merger mechanism for these 
unusual objects. One of the systems, SDSS J105353. 89+520031.0 has a 70% 
chance of having a low-mass white dwarf companion. This system will probably 
form a single helium-enriched subdwarf O star. All four white dwarf systems 
have unusal mass ratios of < 0.2 — 0.8 that may also lead to the formation of AM 
CVn systems. 
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Subject headings: stars: low-mass — white dwarfs — stars: individual (SDSS 
J082212. 57+275307.4, J084910.13+044528.7, J105353.89+520031.0, J143633.29+501026. 



1. INTRODUCTION 

Mergers of binary white dwarfs (WDs) have been proposed to explain supernovae (SNe 
la events, extreme h elium stars including R Coronae Borealis (RCrB) stars, and single sub d 
warf B and O stars Jlben & Tutukovlll984! : IWebbinlJ[l9~8ll : ISaio & Jeffen]l200ol . I2OO2I : iHeber 
20091 ). However, radial velocity surveys of WDs have not revealed a large binary population 



2001. 2002; Karl et a 



2003 



that will merge within a Hubble time ( Marshl Il995t iMaxted et al.l l2000l : iNapiwotzki et al. 



Ne 



emans et al. 



200 51) . In additio n to a few pre- WD + WD 



merger systems (e.g. iGeier et al.l 120071 : iTovmassian et al.ll2010f ). INapiwotzki et al.l ( )2004af ) 
identify only eight merger candidates from the SN la Progenitor Survey (SPY) and the 
literature. 

Radial velocity observations of extremely low-mass (ELM, M ~ 0.2 M Q ) WDs provide 
a new opportunity to find short period binaries. The Universe is not old enough to produce 
ELM WDs through single star evolution. These WDs must therefore undergo significant 
mass loss during their formation in binary systems. The majority of ELM WDs have been 
identified as companions to milli-second pulsars. However, not all ELM WDs have such 



companions (Ivan Leeuwen et al.l 120071 ; lAgueros et al.ll2009al ). Radial velocity, radio, and x- 
ray observations of the lowest gravity WD found in the Sloan Digital Sky Survey (SPSS ), 
J0917+4638, show th at the companion is almost certainly another WD (jKilic et al.ll2007al lb; 
Agiieros et al1l2009bh . 



Recently, lEisenstein et al.l ( 120061 ) identified a dozen ELM WDs in the SDSS Data Re- 
lease 4 area. Here, we present radial velocity observations of four WDs from that sample; 
SDSS J082212.57+275307.4, SDSS J084910. 13+044528.7, SDSS J105353.89+520031.0, and 
SDSS J143633. 29+501026. 8. Our observations are discussed in Section 2; an analysis of the 
spectroscopic data and the nature of the companions are discussed in Section 3. The future 
of these binary systems and the merger products are discussed in Section 4. 



1 Based on observations obtained at the MMT Observatory, a joint facility of the Smithsonian Institution 
and the University of Arizona. 
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OBSERVATIONS 



We used the 6.5m MMT equipped with the Blue Channel Spectrograph to obtain mod- 
erate resolution spectroscopy of four ELM WDs on UT 2009 March 27-29 and April 1-3. In 
addition, we observed J0822+2753 on UT 2008 September 23-24. We used a 1" slit and the 
832 line mm -1 grating in second order to obtain spectra with a wavelength coverage of 3600 
— 4500 A and a resolving power of R = 4300. We obtained all spectra at the parallactic 
angle and acquired He-Ne-Ar comparison lamp exposures after every science exposure. We 
checked the stability of the spectrograph by measuring the centroid of the Hg emission line at 
4358. 34A. The line is stable to within 10 km s _1 , with an average offset from the rest wave- 
length of —0.4 ± 4.9 km s" 1 . We flux-calibrated the spectra using blue spectrophotometric 
standards (IMassey et al.lll988[ ). 



Brown et al.l (120061 ) observed J1053+5200 in 2006 February as part of their hypervelocity 



B-star survey. We include this additional MMT spectrum in our analysis to ext end the time 
baseline. We use the cross- correlation package RVSAO ( IKurtz fc Mink! Il998l ) to measure 
heliocentric radial velocities. We obtain preliminary velocities by cross-correlating the ob- 
servations with bright WD templates of known velocity. However, greater velocity precision 
comes from cross-correlating the objects with themselves. Thus we shift the individual spec- 
tra to rest-frame and sum them together into a high signal-to-noise ratio template spectrum 
for each object. Our final velocities come from cross-correlating the individual observations 
with these templates, and are presented in Table 1. 



Table 1. Radial Velocity Measurements 



Object 


HJD 


Heliocentric Radial Velocity 




+2450000 


(km s- 1 ) 


J0822+2753 


4732.98282 


-254.7 ± 13.6 




4732.99047 


-303.4 ± 11.9 




4733.00223 


-312.4 ± 13.5 




4733.97631 


-295.7 ± 13.4 




4733.99015 


-238.0 ± 14.1 




4917.61018 


+216.2 ± 10.8 




4922.62095 


-291.8 ± 4.5 




4922.68488 


+99.5 ± 4.8 




4922.75956 


+97.1 ± 7.6 




4922.79519 


-121.5 ± 7.8 




4923.62620 


-148.0 ± 5.3 




4923.66669 


+128.6 ± 6.4 




4923.70522 


+223.9 ± 7.3 




4923.76270 


-60.8 ± 5.9 




4924.62160 


-37.1 ± 8.3 




4924.74512 


-99.8 ± 11.2 


J0849+0445 


4918.64083 


+415.3 ± 11.0 




4922.64106 


+247.1 ± 14.1 




4922.70390 


-195.9 ± 10.4 




4922.77816 


-283.1 ± 11.2 




4923.64966 


-144.7 ± 12.9 




4923.73361 


-12.6 ± 10.9 




4923.78104 


-136.6 ± 12.7 




4924.64406 


-34.2 ± 11.0 




4924.66382 


-292.4 ± 13.5 




4924.68169 


+189.2 ± 18.8 




4924.70222 


+391.2 ± 10.3 




4924.72610 


-142.9 ± 14.9 


J1053+5200 


3790.79575 


-45.2 ± 8.5 




4919.62596 


+287.6 ± 8.9 




4922.66788 


-261.0 ± 9.3 




4922.74298 


+98.6 ± 10.4 
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We also use WD model spectra with atmospheric parameters customized for each object 
(see §3) to measure radial velocities. The results are consistent within 10 km s _1 . Thus, 
the systematic errors in our measurements are < 10 km s -1 ; the mean velocity difference 
between the analyses is 1.5 ± 3.7 km s _1 . This small uncertainty gives us confidence that 
the velocities in Table 1 are reliable. 



3. RESULTS 

All four targets display radial velocity variations of > 540 km s _1 between different 
observations. We weight each velocity by it s associated error and solve for the best-fit orbit 



using the code of iKenyon fc Garcia I ( 119861 ). The heliocentric radial velocities are best fit 
with circular orbits and with velocity semi-amplitudes K = 265 — 367 km s -1 . The best- 
fit orbital periods range from 0.0426 to 0.2440 days (1.0 to 5.9 hr). Figures 1-4 show the 
observed radial velocities and the best fit orbits for our targets. We have velocity data from 
3 — 6 nights spread over a time baseline of 3 — 1134 nights. The short orbital periods and 
relatively long time baseline help us to constrain the orbital periods accurately. With orbital 
periods of 1.0 — 1.1 hr, J1053+5200 and J1 436+5010 are the shortest period binary WDs 



yet discovered (see also iMullally et al.ll2009f ). 



The orbital parameters for our targets, including the orbital period, semi-amplitude (K) 
of the radial velocity variations, systemic velocity, the time of spectroscopic conjunction, 
and mass function, are presented in Table 2. The systemic velocities in this table include 
a gravitational redshift term, which should be subtracted from these velocities to find the 
true systemic velocities. This correction is on the order of 3 km s _1 for our targets (see the 
discussion below). 

We perform model fits to each individual spectrum and also to the average composite 
spectra using synthetic DA WD spectra kindly provided by D. Koester. We use the individual 
spectra to obtain a robust estimate of the errors in our analysis. Figure 5 shows the composite 
spectra and our fits using the entire wavelength range. The best-fit T e g- and log g values are 
given in Table 3. We obtain best-fit solutions of 8880 — 16550 K and logg = 6.23 — 6.69 for 
our targets, confirming that they are ELM WDs. 

Slight differences between the continuum level of the observations and that of the best-fit 
model spectrum are evident for J1053+5200. These differences are suggestive of an imperfect 
flux calibration. If we normalize (continuum- correct) the composite spectra and fit just the 
Balmer lines, we obtain best-fit solutions that differ by 70 — 800 K in T e s and 0.1 dex in logg 
for our four targets. These fits are shown in the right panel of Figure 5. These solutions are 
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Table 1 — Continued 



Object HJD Heliocentric Radial Velocity 

+2450000 (km s" 1 ) 



4922.82846 


+10.8 ± 7.4 


4922.91686 


-65.2 ± 7.8 


4923.68895 


-220.1 ± 7.8 


4923.80370 


+177.3 ± 7.7 


4923.87868 


+220.5 ± 6.4 


4924.76454 


-45.3 ± 14.0 


4924.78045 


+260.2 ± 13.7 


4924.79571 


-249.6 ± 14.7 


4924.81573 


+235.9 ± 13.5 


4924.83218 


-60.0 ± 13.5 



J1436+5010 4922.85253 


+301.4 ± 7.4 


4922.93341 


-50.8 ± 11.5 


4922.98659 


+224.2 ± 5.6 


4923.82133 


+230.8 ± 6.8 


4923.85329 


+89.8 ± 6.7 


4923.89751 


+41.9 ± 5.8 


4923.93904 


-188.1 ± 7.0 


4923.98171 


-296.2 ± 6.6 


4923.99118 


+109.6 ± 6.1 


4924.94756 


-140.7 ± 9.3 


4924.95335 


+112.2 ± 11.8 


4924.96121 


+327.2 ± 10.6 


4924.96700 


+243.6 ± 15.9 


4924.97398 


-72.7 ± 11.0 


4924.97975 


-324.1 ± 11.0 


4924.98684 


-372.5 ± 14.0 



Table 2. Orbital Parameters 



Object 


P 


K 


Systemic Velocity 


Spectroscopic Conjunction 


Mass Function 




(days) 


(km s -1 ) 


(km s- 1 ) 


(HJD) 




J0822+2753 


0.2440 ± 0.0002 


271.1 ± 9.0 


-52.2 ± 4.5 


2454732.8312 


0.5038 ± 0.07745 


J0849+0445 


0.0787 ± 0.0001 


366.9 ± 14.7 


47.8 ± 7.4 


2454918.6200 


0.4026 ± 0.06380 


J1053+5200 


0.0426 ± 0.0001 


264.8 ± 15.0 


21.4 ± 7.7 


2453790.7977 


0.08195 ± 0.01418 


J1436+5010 


0.0458 ± 0.0001 


347.4 ± 8.9 


-30.2 ± 5.1 


2454922.8430 


0.1990 ± 0.02969 
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Fig. 1. — The radial velocities of the white dwarf J0822+2753 observed in 2008 September 
(top left panel), 2009 March (middle left panel) and 2009 April (bottom left panel). The 
right panel shows all of these data points phased with the best-fit period. The solid line 
represents the best-fit model for a circular orbit with a radial velocity amplitude of 271.1 
km s -1 and a period of 0.2440 days. 



J0849+0445 




923-7 923.8 924.6 924.7 0.5 1 1.5 

Julian Date + 2454000 Phase 



Fig. 2. — The radial velocities of the white dwarf J0849+0445 observed in 2009 March and 
April (left panels). The right panel shows all of these data points phased with the best-fit 
period. The solid line represents the best-fit model for a circular orbit with K = 366.9 km 
s- 1 and P = 0.0787 days. 
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J1053+5200 




Julian Date + 2453000 Phase 



Fig. 3. — The radial velocities of the white dwarf J1053+5200 observed in 2006 February, 
2009 March and April (left panels). The right panel shows all of these data points phased 
with the best-fit period. The solid line represents the best-fit model for a circular orbit with 
K = 264.8 km s^ 1 and P = 0.0426 days. 




Fig. 4. — The radial velocities of the white dwarf J1436+5010 observed in 2009 April (left 
panels). The right panel shows all of these data points phased with the best-fit period. 
The solid line represents the best-fit model for a circular orbit with K = 347.4 km s _1 and 
P = 0.0458 days. 
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Fig. 5. — Spectral fits (red solid lines) to the composite spectra of our targets (jagged lines, 
left panels) and to the flux- normalized line profiles (right panels). 



Table 3. Physical Parameters 



Object 


T e f{ 


log9 


Mass 


Companion Mass 


i = 60° 


NS 


SN la 


Merger Time 




(K) 




(M e ) 


(Mq) 


(Mq) 


Probability 


Probability 


(Gyr) 


J0822+2753 


8880 ± 60 


6.44 ± 0.11 


0.17 


> 0.76 


1.05 


18% 


5% 


< 8.42 


J0849+0445 


10290 ± 250 


6.23 ± 0.08 


0.17 


> 0.64 


0.88 


15% 


4% 


< 0.47 


J1053+5200 


15180 ± 600 


6.55 ± 0.09 


0.20 


> 0.26 


0.33 


4% 


1% 


< 0.16 


J1436+5010 


16550 ± 260 


6.69 ± 0.07 


0.24 


> 0.46 


0.60 


9% 


4% 


< 0.10 
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consiste nt with the fits to the entire spectra within the errors. Our results are consistent 
with the lEisenstein et al.l (120061 ) analysis within 500 K in temperatur e and .4 dex in surface 
gravity. Similarly, they are also consistent with the iMullally et al.l ( 120091 ) analysis within 
1000 K in temperature and 0.2 dex in surface gravity. Given the higher resolution and 
higher signal-to-noise ratio MMT data, shorter exposure times, and extended blue coverage 
that includes gravity- sensitive higher order Balmer lines, our T e g, log g, and orbital period 
estimates (Table 2 and 3) should be more reliable. 

Figure 6 shows the effective temperatures and surface gravities for our targets (red cir- 
cles) plus the previously identified ELM WDs in the literature. Open c ircles show the WD 



companions to milli-sec ond pulsars PSR J1012+5307 and J1911-5958A (Ivan Kerkwiik et al. 



1996; iBassa et al.ll2006l ). Filled triangles show the ELM sd B star HD 188112 flHeber et al. 



2006; 



2003]) , and the WDs SPSS J0917+4638 fcilic et al.l l2007allbh and LP400-22 fcawka et al' 



Kilic et al.l hood , IVennes et aD I2OO9I ) . All of these WDs show radial velocity 



varia- 



tions. Solid lines show the constant m a ss tra cks for low mass WDs from our updated model 
calculations based on the iPanei et al.l (120071 labeled in M on the right side of the figure) 
study. We mode l the evolution of He -core WDs in clos e bina ry systems, improving on the 
earlier mo dels of lAlthaus et al.l (boOll ) and IPanei et al.l ( 120071 ). Instead of M > O.18M as 



found in lAlthaus et al.l (1200 ll ). we find that for masses M > 0.17017 M diffusion- induced 
hydrogen-shell flashes take place, which yield small hydrogen envelopes. The models with 
M < 0.17017 M Q do not experience thermonuclear flashes. As a result, they have massive 
hydrogen envelopes, larger radii, lower surface gravities, and they evolve much more slowly 
compared to more massive WDs. 

Figure 6 shows that both J0822+2753 and J0849+0445 are 0.17 M WDs, whereas 
J1053+5200 and J1436+5010 are 0.20 M and 0.24 WDs, respect i vely. The errors in 
these mass estimates are ~ 0.01 M . The models by lAlthaus et al.l ( 1200 ll ) also predict 
masses within 0.01-0.02 M of the above estimates. 



4. DISCUSSION 



4.1. J0822+2753 



The temperature and surface gravity for J0822+2753 imply an absolute magnitude 
of M g = 10.1 mag, a radius of 0.04 R & , and an age of 1.2 Gyr. This luminosity places 
it at a distance of 430 pc, 250 pc above the plane. J0822+2753 has a proper motion of 
fi a cos5 = 3.4 ± 3.5 and jis = —19.2 ± 3.5 mas yr -1 (IMunn et al.ll20Q4i ) . Based on the mass 
and radius estimates, the gravitational redshift of the WD is 2.7 km s . Therefore, the true 
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25000 20000 15000 10000 

T eff (K) 



Fig. 6. — The best fit solutions for the surface gravity and temperatures of our targets 
(filled circle s ), ove rlaid on tracks of constant mass from our new calculations based on the 



Panei et al.l (120071 ) models. Spectroscopically confirmed ELM WDs and sdBs found in the 



literature are shown as filled triangles. Companions to milli-second pulsars PSR J1012+5307 
and J1911-5958A are shown as open circles. 
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systemic velocity is —54.9 km 1 . The velocity components with respect to the local standard 
of rest as defined bv lHogg et aD hood ) are U = 64 ± 6, V = - 22 ± 8, and W = -25 ± 7 km 



s _1 . Clearly, J0822+2753 is a disk star f lChiba k Beersl 120001 ) 



We combine the spectra near maximum blue-shifted radial velocity and near minimum 
radial velocity into two composite spectra. If there is a contribution from a companion 
object, it may be visible as an asymmetry in the line profiles. We do not see any obvious 
asymmetries in the line profiles and our optical spectroscopy does not reveal any spectral 
features from a companio n object. A main-seq uence star companion with M > 0.76 M & 
would have Mi < 6.5 mag ( jKroupa fc Toutlll997l ). brighter than the low-mass WD (Mj « 9.3 
mag) and detectable in the /-band. 

Figure 7 shows the SDSS photometry of our targets compared to WD model predictions. 
None of the targets in our sample, including J0822+2753, shows excess flux in the optical. 
Hence, a main-sequence star companion is ruled out for J0822+2753. 

Using the mean inclination angle for a random stellar sample, % = 60°, the companion 
has a mass of 1.05 M & and a separation of 1.8 R & . This separation is about 40 x larger than 
the r adius of the WD. The p robability of a neutron star companion, i.e. M = 1.4 — 3.0 M , is 



18%. lAgiieros et al.l ( l2009af ) obtained radio observations of all four of our targets, including 



J0822+2753, using the Green Bank Telescope. They do not see any evidence of a milli- 
second pulsar companion. Therefore, the companion is most likely a massive WD. However, 
follow-up x-ray observations are required to rule out a neutron star companion. 



4.2. J0849+0445 

J0849+0445 falls in between the cooling sequences for M = 0.17017 and 0.1705 M & in 
the temperature versus surface gravity diagram (Figure 6). Thus, J0849+0445 is 1-2 Gyr 
old, has M g fa 9.5 mag, and is located at 930 pc away from the Sun. The relatively large 
uncertainty in the age estimate is due to the differences between the models with masses 
larger /sm aller than 0.17017 M . The measured proper motion for J0849+0445 is consistent 



with zero (IMunn et al.l 12004] ) . The gravitational redshift of the WD is f»2.3 km s 1 . Hence, 
the true systemic velocity is 45.5 km s _1 . The velocity components with respect to the local 
standard of rest are U = -21 ± 13, V = -24 ± 13, and W = 25 ± 14 km s -1 . J0849+0445 is 
also a disk star. 

As in the case of J0822+2753, we do not see any evidence of a companion in our spectra. 
Main-sequence companions are ruled out based on the SDSS photometry. In addition, radio 
observations do not reveal a milli-second pulsar companion. There is a 70% probability that 
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Fig. 7. — Spectral energy distributions of our targets (error bars) and the WD model pre- 
dictions (open circles). 
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the companion is less massive than 1.4 M®. Therefore, the companion is most likely another 
WD with M > 0.64 M . Assuming an inclination angle of 60°, the companion is probably 
a 0.88 M Q object at 0.8 R® (or 15 x Rwd) separation. 



4.3. J1053+5200 

J1053+5200 has M = 0.20 M®,R = 0.04 R & , and M g = 8.7 mag. This absolute 
magnitude corresponds to a WD cooling age of 160 Myr. Its distance is 1.1 kpc, 1 kpc 
above the plane. It has a proper motion of /j, a cos5 = —29.7 ± 3.5 and fig = —31.2 ± 3.5 mas 



yr 1 ( IMunn et al.l 120041 ) . The gravitational redshift of the WD is 3.2 km s 1 . The velocity 
components with respect to the local standard of rest are U = —110 ± 19, V — —192 ± 28, 
and W = 8 ± 12 km s -1 . J1053+5200 lags by ~ 200 km s" 1 behind the Galactic disk and is 
clearly a halo star. 

The mass function for J1053+5200 implies that the companion is more massive than 
0.26 M . A 0.26 Mq main-sequence star would be about 50% fainter than the WD in the 
/—band and it would have been detected in the SDSS i— and z— band data. The lack of 
excess flux in the SDSS photometry (Figure 7) rules out a main-sequence companion. A 
neutron star companion is unlikely (4% probability). The companion is most likely another 
WD, and specifically a low-mass (M < 0.45 M ) WD. J1053+5200 has a 70% chance of 
having a low- mass WD companion. Assuming an inclination angle of 60°, the companion is 
probably a low mass object with M = 0.33 M & at 0.4 R & separation. 



4.4. J1436+5010 

J1436+5010 is a 40 - 70 Myr old M « 0.24M o WD at a distance of « 800 pc, 710 
pc above the Galactic plane. The gravitational redshift of the WD is 4.2 km s" 1 . It has a 



proper motion of fi a cos5 = 7.8 ± 3.5 and fig = —5.1 ± 3.5 mas yr 1 ( jMunn et al.ll2004j ). The 
velocity components with respect to the local standard of rest are U = 44 ± 14, V = —5 ± 12, 
and W = -28 ± 8 km s -1 . J1436+5010 is a disk object. 

We do not see any evidence of a companion in our MMT spectra of this object. In 
addition, the companion has to be more massive than 0.46 M Q , and such a main-sequence 
star companion is ruled out based on the SDSS photometry data. There is a 9% chance that 
the companion is a neutron star. However, the companion is most likely a carbon/oxygen 
core WD. Assuming an inclination angle of 60°, the companion is probably a 0.60 M & object 
at 0.5 R® separation. 
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Ramirez et al.l ( 120071 ) used Gaussian probability distributions to assign stars to the 



different Galactic components based on kinematics. Using the same criterion, J0822+2753, 
J0849+0445, and J1436+5010 have 94-97% chances of being thin disk members. On the 
other hand, J1053+5200 has a 99% chance of being a halo member. These statistics confirm 
our membership assignments presented above. 



THE COMMON-ENVELOPE PHASE 



Our radial velocity observations and the available optical photometry show that none 
of our targets have main-sequence companions. The probability of neutron star companions 
(M > 1.4 Mp)) ranges from 4% to 18%. However, no such companions are visible in the 
radio data (lAgiieros et al.ll2009al ). Therefore, the companions are most likely other WDs. 



The formation scenarios for close WD pairs include two consecutive common-envelope 
phases or an Alg ol-like stable mass transfer phase followed by a common-envelope phase (e.g. 
Iben et al.lll997l ). The mass ratios (q = Mb r i g ht/-/Wdi m ) for our targets r ange from 0.22 to 0.77 



and they favor a scenario involving two common-envelope phases (see iNelemans et al.ll2000 
and references therein). 



Nelemans et al.l (120001 ) and INelemans &: Toutl (120051 ) demonstrate that the standard 



common-envelope formalism (the a-formalism equating the energy balance in the system) 
does not always work. Instead, they suggest that the common-envelope evolution of close WD 
binaries can be reconstructed with an algorithm (7-algorithm) imposing angular momentum 
balance. Studying the p r ior ey olution of 10 WD+WD binaries where both WD masses are 
known, INelemans et al.l (120051 ) find that the 7-algorithm is able to explain the observed 
properties of those systems following two common-envelope phases. Most of the systems can 
be accounted for with a single-valued 7 = 1.5, where 7 is the rate of angular momentum 
loss as defined by iPaczyhski &: Ziolkowskil (119671 ) . However, 7 is not a stiff parameter, and 
it may vary between 0.5 and 4. The same sy stems may be explained by means of various 
7s and assuming different initial parameters ( INelemans fc Toutl 120051 ). The mass ratios of 
these 10 systems are on the order of unity. 



Kilic et al.l (j2007bl . 120091 ) used the 7-algorithm to study the prior evolution of the ELM 



WD systems J0917+4638 and LP400 -22. These two syst ems have mass ratios of < 0.61 and 
< 0.46, respectively. Using 7 = 1.5, kilic et al.l fcpQ7bh find that the J0917+4638 system 
can be explained as the descendant of a binary system including a 2.2 and a 0.8 M star at 
an orbital separation of 0.4 AU. A simila r study of the LP 400-22 system shows that 7 = 2 
may be more appropriate for this binary (IKilic et al.l 120091 ) . 
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In order to understand the prior history of the short-period binary systems presented in 
this paper, we consider an initial system involving a 3 M and a 1 M star. The 3 M star 
evolves off the main- sequence, overflows its Roche lobe as a giant with a 0.6 M & core, forms 
a helium star (sdB) which does not expand after He-exhaustion in the core, and turns into 
a WD. The 1 M star also overflows its Roche lobe after main-sequence evolution, when its 
core is 0.24 M Q . We also assume that the current orbital period is 0.0458 days (as in the 
J1436+5010 system). The current separation of this system is 0.5 R & . Evolving the system 
back in time, and considering the la test common-envelope phase, the 1 M & star had a 0.24 
M core wh en its radius wa s 10 R Q (llben &: Tutukovlll985[ ). Taking this radius as the Roche 
lobe ra dius (Egglcton 1983 ), th e separation prior to the second common-envelope phase is 
25 R & ( Nelemans &: Tout[2005 ). which corresponds to 7 = 2.01. 

Considering the first common-envelope phase, the 3 M star has a 0.6 M & core when 
its radius is 410 Rq. Taking this radius as the Roche lobe radius, the separation prior to the 
first common-envelope phase is 860 R Q (4.0 AU), which corresponds to 7 = 1.58. Similar to 
this sytem, the common-envelope evolution of all four systems discussed in this paper can 
be explained with 7 = 1.6 — 2.2. In this scenario, the companions are older, smaller, and 
fainter than the 0.17 — O.24M WDs observed today, consistent with the lack of evidence for 
the presence of companions in the SDSS photometry and our optical spectroscopy. 



6. THE FUTURE: MERGER PRODUCTS 

Short period binaries may merge within a Hubble time by losing angular momentum 
through gravitational radiation. The merger time for such binaries is 

(Mi + M 2 ) 1/3 «/, 7 



where the masses are in solar units and the period is in hours (lLandau fc Lifshitzlll958l ). For 
minimum mass companions (i = 90°) the merger times for our targets (in right ascension 
order) are 8420 , 470, 160, and 100 Myr, respectively. All four targets will merge within a 
Hubble time. We now explore possible outcomes of the merger process. 



6.1. EXTREME HELIUM STARS 



Hydrogen- deficient luminous stars, in order of decreasing effective temperature, include 
extreme helium, RCrB, and hydrogen-deficient carbon stars. Studies of the chemical compo- 
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sitions of these stars suggest that they form an evolutionary sequence (I Garcia- Hernandez et al 
20091 ). There are two leading scenarios to explain the origin of extreme helium stars. In 
one scenario, the merger of a He-star o r a He-core WD with a carbon/oxygen core WD 



form s a hydrogen-deficien t supergiant ( jWebbink! Il984l ; llben &: Tutukovl 11984 ; llben et al. 



19961 ; ISaio fc Jefferyl 120021 ). The other scenario, commonly referred to as the born-again 



scenario, suggests that a hydrogen-deficient s tar forms when a post-asymptotic giant branch 
star experiences a late helium shell flash (see llben et al.lll983l . and references therein). This 
flash converts the hydrogen-rich envelope to helium, creating a hydrogen-deficient star. 

Studies of elemental and isotopic abundances for carbon, nitrogen, and oxygen are 
useful for diffe r entiat ing b etween these two sce narios. Based on the observed abundances, 
Saio & Jefferyl (120021 ) and iPandey et all ( 120061 ) argue that most hydrogen- defic i ent ca rbon 
stars and RCrB stars form through WD mergers. In addition, IClayton et all ( 120071 ) find 
that every RCrB and hydrogen-deficient carbon star that they have observed has enhanced 
18 0/ 16 ratios compared to the solar value. They propose the W D merger scen a rio as 
a likely formation mechanism. Based on preliminary calculations, IClayton et al.l (120071 ) 
suggest that the accretion of the helium WD by the carbon/oxygen WD is rapid and it 
induces nucleosynthesis, which converts 14 N to 18 by a-capture. 

The overproduction of 18 O is not predicted in the born-a gain scenario, because either th e 
14 N is burnt to 22 Ne or 18 is destroyed by proton capture (iGarcfa-Hernandez et al.l 120 09). 



There is at least one staj§ identified as a p roduct of the born-again scenario (see IClayton et al.l 
20061 ; iMiller Bertolami fc Althaud 120071 and references therein), Sakurai's object (V4334 
Sgr). Therefore, not all extreme helium stars form through WD mergers. 

Three of the ELM WDs in our sample, J0822+2753, J0849+0445, and J1436+5010, have 
companions more massive than 0.46 M . These 3 systems most likely have carbon/oxygen 
WD companions. They will merge within 100 Myr to 8.4 Gyr and form extreme helium 
stara^l Therefore, these 3 systems provide independent evidence that there is a mechanism 
to form hydrogen deficient stars through WD mergers. 



6.2. SINGLE HOT SUBDWARF STARS 

Close binary evolution plays an important role in the formation of subdwarf B stars, 
as witnessed by the large fraction (> 40%) of sdB stars in binaries. Mass transfer between 



2 Two more have been discussed extensively in the literature, FG Sge and V605 Aql. 

3 With a mass ratio of < 0.22, J0822+2753 may become an AM CVn system instead (see Section 6.5). 
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the companions and common-envelope evolution can lead to large a mount of mass loss prior 
to the start of core He-burning, creating an sdB star in a binary. ICastellani fc Castellani 
( 119931 ) demonstrated that if low-mass stars lose enough mass on the red giant branch, they 
can depart the red giant branch before the core He flash. Instead, the He flash happens on 
the hot WD cooling track. Depending on when this flash happens, sdB or He-enriched sdO 
stars form. 



Only about 4% of He-enriched sdO stars are in binary systems (INapiwotzki et al.ll2004bl ). 
The mergers of two He-core WDs can explai n the lower fraction of binaries observed among 
the He-enriched sdO stars (IHan et all 120031 ). With shrinking s eparation, th e less massive 
object is accreted onto the companion, leading to He ignition (jHeberl 120091 ) . The merger 
produ ct will be enriched i n CNO, similar to the observed abundances of He-enriched sdO 
stars JSaio fc Jefferviboooh . 



In addition to the merger scenario, single hot subdwarf stars can also be produced 

through a common-envelope phase with a massive planet or a brown dwa rf. The recent 

discoveri es of planet s arou nd the subdwarf B stars V Pegasi (ISilvotti et all 120071 ) and HD 
149382 (jGeier et al.l 120091 ) show that this channel of formation contributes to the single 
subdwarf B star population. 

There is a 70% chance that the companion to one of our targets, J1053+5200, is a 
low- mass He-core WD. The merger of this binary WD system will most likely create a He- 
enriched sdO star. Hence, binary mergers of two He-core WDs contribute to the single hot 
subdwarf population in the Galaxy. 



6.3. SNe la and .la 



SNe la are caused by the thermonuclear explosion of WDs growing to the Chandrasekhar 
mass either by accretion from a companion or by mergers of two degenerate stars. The 
double degenerate scenario requires mai nly mergers of two CO WDs to have su fficient mass 
to exceed the Chandrasekhar mass limit (jWebbinklll984! ; Ilben fc Tutukovlll984l ). However, if 
the mass transfer is stable, accretion from a He-core WD can also result in accumulation of 
Chandrasekhar mass by the C O WD accretor a nd this can result in a SN la explosion. The 
population synthesis models by lYungelsonl ( 120051 ) demonstrate that the expected contribution 
of He+CO WD binaries to the SNe la rate is two orders of magnitude smaller than that of 
the CO+CO WD systems. 



The visible components of the four binaries discussed in this paper are 0.17-0.24 M & 
WDs. There is a 1-5% chance that the companions are massive WDs and the total mass 
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of the binary systems exceed 1.4 M & (see Table 3). Our targets are most likely not SNe la 
progenitors. 



Recently, iGuillochon et al.l ( 120091 ) presented a new mechanism for the detonation of a 
sub-Chandrasekhar mass CO WD through accretion from a low- mass He-core WD. If the 
mass accretion is dynamically unstable, the instabilities in the accretion stream can lead 
to the detonation of surface helium that accumulates on the CO primary during the final 
few orb its prior to merger. These detonations are likely to resemble dim Type la SNe 
(or .la, iBildsten et al.l 120071 ) . and would primarily synthesize intermediate- mass elements. 
Under certain conditions, the primary itself is also detonated. An important feature of this 
mechanism is that the total system mass does not need to exceed the Chandrasekhar limit. 
As such, the systems we present in this paper could very well be the progenitor systems for 
these events. 



6.4. WHITE DWARF + NEUTRON STAR MERGERS 



Green Bank Telescope observations do not detect milli-second pulsar companions around 
our targets. However, neutron star companions cannot be ruled out based on the radio 
data alone. X-ray observations can detect the blackbody emission from a neutron star 
companion even if it is radio-quiet or if its pulsar beam misses our line of sight. XMM- 
Newton observations of the extremely low mass WD, SDSS J0917+4638, do not detect a 
neutron star in the system (lAgiieros et al.ll2009bl ). Similar observations will be necessary to 
search for neutron star companions in our targets. 



Nelemans et al.l (120011 ) estimate that NS + WD merger rate is two orders of magnitude 



smaller than WD + WD merger rate. However, based on the mass fun ction alone, there is 



a 4-18 % chance that the companions to our targets are neutron stars. iGarcia-Berro et al. 



(120071 ) study the evolution of a 1.4 M & neutron star and a 0.6 M & merger using a smooth 
particle hydrodynamics code. They find that, once the white dwarf has filled its Roche 
lobe, it is disrupted in a few orbital periods, e.g. on the order of 5 minutes. The final 
configuration consists of a neutron star surrounded by an accretion disk, and the mass loss 
from the system is negligible. Ho wever, the outcome of mergers of NS with lighter WDs is 

(2002) suggest that mass exchange in a 0.2 + 1.4 M & system may 
(120 lOl ) argue for WD or He-star donors in several ultra-compact 
x-ray binary systems. These donors have initial masses larger than 0.32-0.35 M & and they 
somehow survive disruption in the merger with the NS. Therefore, 0.2 M WDs may also 
survive mergers with neutron stars. 



not clear. lYungelson et al. 



be stable. Nelemans et al. 
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6.5. AN ALTERNATIVE FUTURE: AM CVn SYSTEMS 



AM CVn stars are interacting double stars with WD accretors and orbital periods less 
than about one hour. There are three formation scenarios for AM CVn systems involv- 
ing three types of donor s tars; WDs, helium stars, or evolved main-sequence stars (see 
Postnov &: Yungelsonl 120061 . a nd references therin). Studying the CNO and He abundances 
of known AM CVn systems, iNelemans et al.l ( 120101 ) find evidence of WD donors in some 
systems, and evolved helium star donors in others. 

The WD channel requires a binary system with short enough orbital period that gravi- 
tational wave radiation drives the stars into contact. The low- mass WD fills its Roche lobe 
and transfers mass to the companion. Depending on the mass ratio of the binary system (if 
the mass ratio is small) , the mass transfer is s table and the system evol yes to longer periods 



(see iMarsh et al.l 12004 ; INelemans et al.l |2010| . and references therein). iMarsh et al.l (12004) 
argue that despite the absence of a single system of extreme mass ratio amongst the observed 
close double WD population, this channel of formation is probable. 

Here, we have uncovered four binary systems with extreme mass ratios of < 0.22 to < 
0.77. The importa nt question is whethe r these systems will merge or if they will instead create 



AM CVn systems. IMarsh et al.l (120041 ) suggest that the mass transfer between double WDs 



Motl et al. 


(2007 


) and 


Racine et al. 


(2007) 



demonstrate that the spin/orbit coupling is strong, raising the critical mass ratio to avoid 
merger from around 0.2 (if there is no coupling) to 0.4-0.7. Based on these studies, we 
would expect the majority of the WDs discussed in this paper to form AM CVn systems. 
However, the SPH and grid-based calculations of mergers of WDs do not completely agree 
on the outcome of c ontact, and the prior evolution of AM CVn systems is still uncertain (see 



Fryer fc Diehlll2008l . for a discussion). Nevertheless, with mass ratios of < 0.27, J0822+2753 



and J0849+0445 are strong candidates for future AM CVn systems. 



7. CONCLUSIONS 



Almost all k nown doubl e WD systems have mass ratios on the order of unity (INelemans et al 



20051 ) . Recently, iKilic et al.l ( 12007bl . 120091 ) identified two WD binary systems with mass ra- 
tios < 0.6. However, these two systems (SDSS J0917 and LP400— 22) will not merge within 
a Hubble time. 



The four binary systems discussed in this paper will merge within a Hubble time. They 
have extreme mass ratios and are likely progenitors of RCrB stars, single He-enriched sdO 
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stars, or AM CVn stars. These systems may even contribute to the SNe la population, but 
the probability of this event is small. 



Liebert et all ( 120051 ) estimate that low-mass (M < O.45M ) WDs make up about 10% 



of the local WD popu lation, corresponding to a formation rate of 4 xlO 14 pc 3 yr 



Eisenstein et all (120061 ) discovered only 13 ELM WDs out of a sample of 9316 WDs found in 
a survey volume of > 4 kpc 3 in the SDSS Data Release 4 footprint. ELM WDs are rare, they 
make up about 0.14% of the local WD sample. Of course, the SDSS spectroscopic sample 
suffers from selection bias and incompleteness, but taken at face value, 0.14% corresponds 
to a formation rate of 6 xl0~ 16 pc~ 3 yr -1 . However, all four ELM WDs discussed in this 
paper are going to merge within a Hubble time, and some within the next few hundred Myr. 
Therefore, the formation rate of ELM WDs may be higher by an order of magnitude. For 
the SDSS DR4 survey volume of 4 kpc 3 , the formation rate is 0.2 — 2 x 10~ 5 yr -1 , with 
the caveat that this n umber is highly uncerta i n due to selection bias present in the SDSS 
spectroscopic survey. iBogomazov fc Tutukovi (120091 ) estimate a CO+He WD merger rate 



of 6-8 xlO 3 yr 1 . The ELM WD mergers discussed in this pa per contrib u te on ly a small 



fraction to the overall CO+He WD merger rate in the Galaxy. lYungelsonl (120051 ) estimate 
a SNe la occurence rate of 10 -5 yr -1 from CO+He WD progenitors. Our estimate of the 
formation rate of ELM WDs is similar to this result. 

In this paper, we present radial velocity data for four stars. Ho wever, we are ob t aining 
radial velocity observations for the remaining ELM WDs in the lEisenstein et all (120061 ) 
sample. To date, almost all of these ELM WDs show significant radial velocity variations 
indicating the presence of a companion star. We are continuing to follow-up these objects at 
the MMT to constrain their orbital periods accurately. The importance of these observations 
is that, we not only find short period binaries that will merge within a Hubble time, but 
also the majority of the binaries have extreme mass ratios. Understanding the prior history 
of these systems requires an understanding of the common envelope phase. Studying the 
mass ratios and physical characteristics of these systems will help in understanding common 
envelope evolution of close binary pairs. 
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